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Recent results for the spectrum of light hadrons provide clear evidence for the failure of quenched QCD and 
encouraging signs that simulations with dynamical sea quarks rectify some of the discrepancies, although string 
breaking has not yet been observed. The use of perturbation theory to match lattice quark masses to continuum 
schemes remains questionable, but non-perturbative methods are poised to remove this uncertainty. The inclusion 
of dynamical sea quarks substantially reduces estimates of the light quark masses. New results for the lightest 
glueball and the lightest exotic hybrid state provide useful input to phenomenology, but still have limited or no 
treatment of mixing. The 0(a)-improved Wilson quark action is well-established in quenched QCD for f3 > 5.7, 
with most parameters obtainable non-perturbatively, in which range scaling violations are small. Progress has 
also been made with high-order improvement schemes for both Wilson and staggered quarks. 



1. INTRODUCTION 

Calculation of the light hadron spectrum is 
an essential part of any numerical simulation of 
QCD, because it provides the most direct way of 
fixing the quark masses. The bare mass for each 
quark flavour is tuned until the masses of corre- 
spondingly flavoured hadrons agree with experi- 
ment, having chosen one dimensionful quantity to 
set the overall scale. It is expensive to simulate 
at masses much below that of the strange quark, 
so the u and d quark masses must be obtained by 
chiral extrapolation. Electromagnetic corrections 
are ignored, so we cannot distinguish the u and d 
masses. 

The hadron spectrum is widely used to test 
scaling for improved lattice actions; the improve- 
ment programme being our best hope for reliable 
simulations with dynamical sea quarks. Once the 
known hadron masses have been used to validate 
the simulations, the spectrum calculations pro- 
vide important predictions for phenomenology. 
They are guiding searches for glueballs and for 
exotic mesons. Finally, the light quark masses 
are some of the poorest-known Standard Model 
parameters and these uncertainties matter; for 
instance, Standard Model predictions of e'/e are 
very sensitive to m s + The challenge 

is to find a reliable way of matching the lattice 



quark masses to those in a continuum perturba- 
tive scheme, which can be used for phenomenol- 
ogy 

1.1. The determination of quark masses 

Quarks do not appear as asymptotic states in 
QCD, because of confinement. As a consequence, 
it is usual to quote values for the running quark 
mass at a particular scale in a specific scheme, eg 
m MS (2 GeV), although, alternatively, we could 
run the mass up to high scales and quote the RG 
invariant (scheme independent) value 
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Given M , it is straightforward to run the mass 
down to any desired scale, using continuum per- 
turbation theory, provided we remain within the 
perturbative regime. 

For Wilson quarks, the pseudoscalar meson 
mass vanishes at K = K C rit and the bare quark 
mass is 
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where n q is determined by tuning the mass of 
a (/-flavoured hadron to its experimental value. 
The renormalised quark mass may be determined 
either using the vector Ward identity (conserved 
current), 
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TO f S (M) = Z s (fta) 1 m q = Z m (pa) m qi 
or the axial Ward identity, 
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which gives the renormalised quark masses with- 
out using K cr it- Here, quantities denoted at the 
scale fj, (a) are renormalised (bare). For stag- 
gered quarks, the remnant chiral symmetry en- 
sures Za = 1, Zp = Zs — Z" 1 , so both methods 
are essentially the same (unfortunately, the 1-loop 
term in the matching is large, casting doubt on 
the reliability of perturbation theory). 

1.2. Some important questions 

In reviewing recent results for the light hadron 
spectrum, I will focus on the following questions. 
Can we conclude that quenched QCD is wrong? 
To what extent have we observed effects which 
can be attributed to dynamical sea quarks? How 
effective are improved actions at reducing scaling 
violations? Is the matching of lattice and con- 
tinuum quark masses under control? Finally, I 
will report on a few phcnomcnologically impor- 
tant predictions for glueballs and hybrids. 

2. QUENCHED APPROXIMATION 

It has been known for some time that the 
quenched approximation gets the spectrum in the 
strange sector wrong and underestimates hyper- 
fine splittings Q . There are several related symp- 
toms of this failure, in particular, using Mk and 
as input gives different values for the strange 
quark mass. 

2.1. Chiral extrapolations 

Chiral extrapolation is now the main source 
of systematic error and a key issue is whether 



the data requires the inclusion of quenched chi- 
ral logarithms. If so, we need to be able to con- 
trol the extrapolation to within a few percent 
to show that the quenched approximation breaks 
down in the ud sector too. CP-PACS has reported 
such results for the Wilson quark action at this 
conference, having increased the statistics at its 
smallest lattice spacing during the year (|^]. As 
yet, though, this failure of the quenched approx- 
imation has not been confirmed using staggered 
quarks, or other improved actions. 

In quenched QCD, the rf meson remains light 
in the chiral limit, because the infinite series of 
loop diagrams, which gives it a non-zero mass 
in the full theory, is absent jjjj. The singlet 2- 
point function has a double pole, giving rise to so- 
called quenched chiral logarithms, which diverge 
as vfiq — > 0. The renormalisation of quenched chi- 
ral perturbation theory has been carried out at 
one loop H . Anomalous chiral behaviour persists 
in dynamical simulations when the valence-quark 
and sea-quark masses are different, as is necessar- 
ily the case for the strange quark in simulations 
with Nf = 2 sea quarks (partial quenching), and 
this can produce strong dependence on the sea- 
quark mass when both it and the valence-quark 
mass are small Q. 

The CP-PACS data are consistent with the 
presence of quenched chiral logarithms 
Both the pseudoscalar meson masses and decay 
constants indicate a non-zero value for the coef- 
ficient of the leading quenched chiral logarithm, 
6 = 0.10(2). Also, using quenched chiral per- 
turbation theory to extrapolate the pseudoscalar 
meson mass to zero, gives a value of K cr it in bet- 
ter agreement with a linear extrapolation of the 
axial Ward identity quark mass, than a simple 
quadratic polynomial in the quark mass. 

2.2. The quenched QCD spectrum 

The evidence of a non-zero value for 5 justi- 
fies CP-PACS's use of quenched chiral perturba- 
tion theory to produce their final results for the 
quenched spectrum. They obtain the continuum 
limit, shown in Fig |Tj, by linearly extrapolating 
data at four lattice spacings, with a -1 ranging 
from 2 to 4 GeV, on lattices which are all roughly 
3 fm in size. 
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Light Hadron Spectrum in Quenched QCD 
final results from CP-PACS: qChPT chiral extrapolations 
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Figure 1. CP-PACS results for the light hadron 
spectrum in quenched QCD obtained using the 
Wilson quark action Mm- 



Fig [l] confirms the picture that there is no 
choice for the strange quark mass which can ex- 
plain the whole spectrum. Using as input, 
the meson hyperfine splitting, the octet baryon 
masses, and the decuplet baryon mass splittings, 
are all too small. Using to fix the strange 
quark mass reduces the discrepancies amongst 
the baryon masses, but the other faults remain. 
The problem is not due to the way in which the 
strange quark mass is defined; as can be seen in 
Fig |2| the vector and axial Ward identity defini- 
tions agree in the continuum limit, but with lim- 
iting values which depend on the choice of input: 

143(6) MeV (Af input) 
115(2) MeV (M K input) 

The CP-PACS result for the u and d quark mass 
is 

mjjjf (2 GeV) = 4.55(18) MeV (9) 
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Figure 2. CP-PACS results for the strange quark 
mass in quenched QCD, defined using the vector 
(Pert in the figure) and axial (WI) Ward identi- 
ties ffl. 



again with consistency between the two defini- 
tions in the continuum limit. The problem with 
the strange quark mass is directly related to the 
low value of J = Mk* (dmv/dm PS ), since 
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if linear chiral behaviour is assumed for m 
m v |. 

It is desirable to confirm the CP-PACS results 
using an improved action. The most direct com- 
parison of this sort is with MILC results |) 10 
for staggered quarks, which are also available on 
lattices of at least 3 fm in size. Chiral extrap- 
olation again turns out to be the most delicate 
issue. For staggered quarks, it is the mass of the 
non-Goldstone pion which appears in quenched 
chiral perturbation theory, and MILC includes 
terms linear in this mass in the chiral extrapola- 
tion of the nucleon and vector meson masses. The 
continuum extrapolation is shown in Fig |^ and, 
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Figure 3. Continuum extrapolation of MILC data 
for m^/m p in quenched QCD, using staggered 
quarks |]|LJ. 

combining various chiral and continuum fits, they 
quote the final result 

— = 1.254 ±0.018 ± 0.027, (11) 

m p 

in agreement with the experimental value of 1.22. 
The statistical error is similar to that of CP- 
PACS, and the difference between the results is 
only 2.5(7, which is not a serious disagreement. It 
would be interesting to combine the MILC data 
with that of Kim and Ohta [|llj , which is on a sim- 
ilar size lattice, but with a smaller lattice spacing 
(f3 = 6.5), to see if this error can be reduced. 

In summary, CP-PACS has provided us with 
evidence of the breakdown of the quenched ap- 
proximation for light hadrons at the few percent 
level, the most striking symptom of this being 
non-uniqueness of the strange quark mass, al- 
though confirmation from improved actions is still 
awaited. 

3. SEA-QUARK EFFECTS 

During the last few years, several groups have 
begun systematically to explore the parameter 



space for simulations with dynamical quarks. 
Finite-size effects must be checked, and this is 
not as straightforward as in quenched QCD, be- 
cause the lattice spacing at fixed sea-quark mass 
depends on the physical quantity used to define 
it, and decreases significantly with the sea-quark 
mass. So the volume is only well-defined for chiral 
sea quarks. 

I will focus on simulations with two degener- 
ate flavours of sea quark, which may represent u 
and d, in which case the strange quark must be 
treated in the quenched approximation. The data 
for Wilson quarks, obtained with unimproved and 
various choices of improved action (SES AM , 
CP-PACS HU and UKQCD Q), correspond 
to mass ratios mps/my between 0.6 and 0.8, and 
lattice sizes of around 2 fm, whereas the MILC 
data for staggered quarks Q extend to a mass 
ratio of around 0.3 on lattices which are some- 
what larger. 

3.1. Wilson quarks 

UKQCD has explored finite-size effects for 
hadron masses, in which valence and sea-quark 
masses are the same, and finds that a significant 
effect occurs at the lower end of the above range 
of quark masses, between lattices sizes of 0.9 and 
1.4 fm [jl4|. Above 1.4 fm, no effect is observed. 
Consequently, chiral extrapolations using linear 
and quadratic polynomials in the valence and sea- 
quark masses, and including cross-terms, is prob- 
ably safe for the Wilson quark data mentioned 
above. 

SESAM notes that data at different sea-quark 
masses are uncorrelated, so that the error in the 
chiral extrapolation is larger than for quenched 
QCD, probably obscuring sea-quark effects 
As reported by CP-PACS the trend in the 

light hadron masses, as lattice spacing decreases, 
is to agree with experiment, but the statistical 
errors are still large compared to the quenched 
QCD results, where the discrepancy with exper- 
iment is small anyway, and CP-PACS does not 
attempt a continuum extrapolation. 

3.2. Staggered quarks 

The results obtained by MILC for N f = 2 stag- 
gered quarks [ [Tcfl are perhaps the most surpris- 
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Figure 4. Edinburgh plot comparing continuum 
results for quenched and Nj — 2 QCD using stag- 
gered quarks, obtained by MILC fief] . 
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Figure 5. Comparison of the strange meson split- 
ting in quenched and Nf — 2 QCD, obtained by 
CP-PACS IBIS. 



ing at this year's conference. They perform a 
const + 0(a 2 ) continuum extrapolation of data at 
five lattice spacings, at fixed m T /m p , to produce 
the Edinburgh plot shown in Fig || Although the 
continuum and chiral extrapolations appear to be 
well behaved, the dynamical results lie systemat- 
ically above the quenched results, and so deviate 
further from experiment. This is true even at 
large quark masses, where the entire calculation 
should be under control! 

3.3. The strange-quark sector 

Returning to the Wilson data, we might hope 
to see an improvement in the strange sector from 
dynamical quarks. CP-PACS reports a tendency 
for the strange meson splitting to increase as the 
continuum limit is approached PJl^|, in better 
agreement with experiment than the quenched 
data, as shown in Fig ||. 

SESAM points out that J no longer has a sim- 
ple interpretation as the slope of the vector me- 
son mass with respect to the square of the pseu- 
doscalar meson mass, mps, because J is no longer 



a function only of mps, but depends separately 
on the valence and sea-quark masses [Q. Thus, 
chiral extrapolation cannot be avoided in comput- 
ing J. At fixed sea-quark mass and fixed lattice 
spacing, there is no improvement in the value of 
J compared to the quenched approximation, and 
there is no discernable dependence on the sea- 
quark mass, as noted also by UKQCD @0. 
This could be due to the continued use of the 
quenched approximation for the strange quark. 

3.4. Quark masses 

Partial quenching, in which the valence and 
sea-quark masses are taken to be different, and 
which is necessary for the strange quark in Nf = 2 
simulations, allows for an alternative definition 
of quark mass at non-zero lattice spacing. The 
partially-quenched quark mass is defined as 

m PQ a =f_i 1 \ (i2) 

where K C rit(K SC a) is defined by the vanishing of 
the pseudoscalar meson mass, at fixed K SC a- Last 
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Figure 6. CP-PACS results for the ud quark mass 
in two-flavour QCD at three lattice spacings, us- 
ing the vector WI (Perturbative in the figure), 
axial WI (WI) and partially-quenched (Perturba- 
tive PQ) definitions |Q. 



year, SESAM reported that K C rit(^sea) depends 
strongly on n sca |]| . 

CP-PACS takes K soa = n u d to define a 
partially-quenched mass. Their results |||l3| in 
Fig || show that, at fixed lattice spacing, the 
partially-quenched, vector and axial Ward iden- 
tity masses all disagree. However, it appears 
that, within the relatively large statistical er- 
rors, all definitions are converging to a contin- 
uum limit of around 2.5 MeV, although the data 
is not good enough to justify an extrapolation. 
As shown in Fig [| CP-PACS also finds that the 
estimates for the strange quark mass, obtained 
using the K and <fi masses, disagree at non-zero 
lattice spacing, but, within the large statistical 
errors, appear to be consistent with a single con- 
tinuum limit around 80 MeV. Evidently, there is 
a large quenching effect which overestimates the 
light quark masses by as much as 40%. It should 
be noted that the results for the strange quark 
mass may still be an overestimate due to the use 
of partial quenching. In fact, CP-PACS estimates 
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Figure 7. Comparison of CP-PACS results for the 
s quark mass, determined from the <f> and from 
the K meson masses, in two-flavour QCD P,p"3[. 



for the ratio m s /m u d from the axial Ward iden- 
tity are essentially independent of lattice spac- 
ing at a value of about 25 |^,[l3|], consistent with 
lowest-order chiral perturbation theory, although 
a slightly smaller value is favoured at next or- 
der |lq| . 

3.5. String breaking 

The much-hoped-for signal of sea-quark effects 
is, of course, string-breaking. Although Philipsen 
and Wittig have demonstrated very clearly the 
crossover from a string-like to a two-meson state 
in a (2 + l)-dimensional SU(2) Higgs model 0, 
zero-temperature QCD results, such as those 
from UKQCD in Fig ||, show little dependence 
of the potential at large distances on sea-quark 
mass |l4j. It is likely that the Wilson loop op- 
erator does not project well onto broken string 
states outside a narrow mixing region ]l7| , and a 
computation of the full matrix correlator of Wil- 
son loops and two-meson operators is needed, but 
see for further discussion of this. 
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4. IMPROVED ACTIONS 

Improvement at least to 0(a) is necessary for 
dynamical quark simulations and considerable 
progress has been made in recent years, partic- 
ularly by the Alpha Collaboration, in implement- 
ing this non-perturbatively for Wilson quarks. 
Most of the scaling tests have been performed in 
quenched QCD, but some dynamical results are 
now becoming available. Higher-order improve- 
ment has the potential for big pay-offs, but many 
parameters have to be fixed reliably. 

4.1. O(a)-improved Wilson quarks 

This requires one additional term in the action, 

a 5 

csw^;^2q(x)i(J^F^(x)q(x), (13) 

X 

as noted first by Shcikholcslami and Wohlert 
(SW), plus explicit determination of the mass de- 
pendence (coefficients bA, ■ ■ • ) and mixing (coef- 
ficients CAi ■ ■ ■) at 0(a) of composite fields: 



C R = Z (l + b am q ){0 + ^ c n aO Tl 



(14) 



Using Schrodinger functional methods at m q = 0, 
Alpha has determined many of the coefficients for 
quenched QCD for (3 > 6.0 (a < 0.1 fm), by im- 
posing chiral symmetry, which is broken at 0(a) 



by the Wilson action (for a review see |19|l ). De 
Divitiis and Pctronzio use the quark-mass depen- 
dence of the PCAC relation to extract bA — bp, 
by - 6s, b m , Z m Z P /Z A and Z m Z s /Z v |20| . 

The difficulty in extending this to larger lat- 
tice spacings is the occurrence of exceptional con- 
figurations. These can be regulated by taking 
m q ^ 0, and the SCRI group finds that the 
mass dependence of csw is very weak, so that 
for (3 > 5.7 ED], 



csw 



1 - 0.6084 g 2 - 0.2015. g 4 + 0.03075.g e 
1 - 0.8743 5 2 



(15) 



within 1% of the Alpha curve for (3 > 6.0. At 
fixed quark mass, defined by my/mpg = 0.7, at 
(3 = 5.7, scaling violations are reduced from 41% 
to 3% in my, and from 33% to 2% in ton 21 1. 

A comparison of the scaling violations in the 
vector meson mass, at the same fixed quark mass, 
for various 0(a)-improved actions is given in 
Fig |^. Consistent continuum extrapolations, lin- 
ear in a 2 , are possible for the non-perturbatively 
improved results, in agreement with 0(a) extrap- 
olation of the Wilson data, but, if the tadpole- 
improved results are also to be consistent, their 
continuum extrapolation must include a term lin- 
ear in a p^ . Also, it is interesting that the stag- 
gered results show big 0(a 2 ) scaling violations. 
The QCDSF Collaboration now has results for 
the non-perturbatively 0(a)-improved Wilson ac- 
tion at a smaller lattice spacing ((3 = 6.4) than in 
Fig |9|, which confirm the smooth approach to the 
continuum limit [25( ]. 

The Alpha Collaboration has computed the 
clover coefficient for non-perturbative 0(a)- 
improvement of two-flavour QCD with Wilson 
quarks, obtaining 



csw = 

1 - 0.454 g 2 - 0.175 g A + 0.012 g 6 + 0.045 g s 
1 -0.720 g 2 

valid for j3 > 5.4 and probably for f3 > 5.2 |2 



(16) 



4.2. Improved staggered quark actions 

The motivation for improving staggered quark 
actions is rather different from that for Wilson 
quarks. The main objective is to reduce flavour 
symmetry breaking by reducing the coupling of 
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ferent 0(a)-improved actions [ggfl . 



quarks to high-momentum gluons. These can 
couple quarks at opposite corners of the Bril- 
louin zone, which correspond to different flavours 
in the continuum. Unfortunately, there are 15 
dimension-6 terms which can contribute to the 
action at 0(a 2 ), making a non-perturbative esti- 
mate of their coefficients almost impossible ]25| . 
One-loop perturbation theory is typically unreli- 
able for staggered quarks, and so the only practi- 
cal option is to use tadpole estimates. 

M1LC has investigated the Naik quark action, 
in which a three- link hopping term is used to can- 
cel 0(a 2 ) terms in the staggered-quark action at 
tree level, and then the coefficients are tadpole- 
improved ^6|. The gluon action they use is 0(a 2 ) 
one-loop tadpole improved. Results are compared 
with those from staggered quarks using both the 
standard, and the improved gluon action. The 
observed scaling violations are 0(a 2 ) and are min- 
imised primarily by the use of improved glue at 
low quark mass, and by the use of the Naik action 



at large quark mass. The speed of light from the 
continuum dispersion relation is closer to unity 
for the Naik action than for the staggered ac- 
tion (both with improved glue), whereas flavour 
breaking, measured from the splitting between 
the non-Goldstone and the Goldstone pions, is 
most sensitive to the use of improved glue, in line 
with the above intuition and previous results us- 
ing 'fat' (ie APE smeared) links (27|. 

Lagae and Sinclair have generalised the M1LC 
fat-link action, achieving a similar reduction in 
flavour violation p8j. For dynamical staggered 
quarks, MILC also finds that improving the gauge 
action or fattening the links improves the flavour 
symmetry [ p9[ . 

Fat links can be used with the SW action to 
allow csw t° be tuned to minimise the spread 
of near-zero modes j30|. The smeared links al- 
low a small enough value of csw to avoid excep- 
tional configurations and, being less sensitive to 
UV modes of the gauge fields, renormalisation 
constants are closer to 1 than for the standard 
SW action. 

4.3. Higher-order improvement 

The most adventurous improved action hunters 
seek small scaling violations on very coarse lat- 
tices. Inspired by the fixed-point action, De- 
Grand has been testing fermion actions which 
couple all the fields on a 3 4 hypercube pj| . He 
finds rotational invariance (ie the speed of light 
from the continuum dispersion relation is close 
to 1) even at a = 0.36 fm, although more com- 
plicated Pauli interactions seem no better than 
the simple clover term, Eq (p"3|), for boosting the 
hyperfine splitting. 

The Q(q 3 )-accurate D234c quark action of Al- 
ford et al. |p2| , with plaquette plus 2x1 rectan- 
gle gluon action and mean-link (Landau gauge) 
tadpole-improved coefficients, has scaling viola- 
tions of only 7% (in at mps/my = 0.7), and 
a much better dispersion relation than the SW 
action, at a lattice spacing as large as 0.4 fm. 

Thus, while non-perturbative O(a) improve- 
ment for Wilson quarks is well-established for 
quenched QCD, and so sets almost a mandatory 
minimum improvement level for dynamical simu- 
lations, there are further encouraging signs that 
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the big gains from higher-order improved actions 
could become a practical reality. 

5. NON-PERTURBATIVE QUARK 
MASS RENORMALISATION 

Most of the existing results for quark masses 
have used perturbative matching and continuum 
extrapolation. There are some indications that 
this procedure may not be fully under control. 
Inconsistencies between the results from using dif- 
ferent definitions and different actions cast doubt 
on the use of perturbation theory, but the issue 
may be clouded by discretisation errors. Conse- 
quently, the main progress this year has not been 
to increase the precision of the mass estimates 
compared to last year's review |33| ], but rather 
in implementing two non-perturbative matching 
schemes. The work described below is entirely in 
the quenched approximation. 

5.1. Perturbative matching 

The CP-PACS results, given in Eqs (§) and 
(||), obtained with the Wilson action and 1-loop 
Z factors, may be compared with the QCDSF 
group's new results using the non-perturbative 
0(a)-improved Wilson action [p3| . 

QCDSF has new data at = 6.2 (on 24 3 x 48 
and 32 3 x 64 lattices) and at j3 = 6.4 (on a 32 3 x 48 
lattice). They determine the quark masses using 
the axial Ward identity, with non-perturbative 
values of Za from the Alpha Collaboration, and 
tadpole-improved 1-loop perturbative values for 
Zp in Eq (0). They fix the strange quark mass 
from Mk and extrapolate results at /3 = 6.0, 6.2 
and 6.4 linearly in a 2 to the continuum, obtaining 



m^f (2 GeV) 

MS 



4.13 ±0.08 MeV 



m, b (2GeV) = 98.1 ±2.4 MeV 



(17) 
(18) 
and 



which are lower than the CP-PACS values, 
significantly so for the strange quark. 

Working at fixed lattice spacing, with tree-level 
0(a)-improved Wilson quarks at the same (3 val- 
ues as QCDSF, Gimenez et al. [p 34 also find 



that the strange quark mass defined by the axial 
Ward identity with perturbative matching is sig- 
nificantly lower than that from the vector Ward 
identity, or from the non-perturbative RI scheme 



(see below). However, the difference could partly 
be due to discretisation effects. 

Before discussing non-perturbative matching, 
it is interesting to note that the first results for 
quark masses using domain wall fermions were 
reported this year Pqj . The preliminary results, 
using Mk to fix the strange quark mass, 
to set the scale, and perturbative matching, at 
(3 = 5.85, 6.0 and 6.3, appear to scale within 
relatively large errors. The weighted average is 
mf S (2 GeV) = 82(15) MeV. While the error is 
too large to add much new information at the mo- 
ment, the result shows that domain wall fermions 
can provide a valuable independent determina- 
tion of quark masses. 

5.2. The non-perturbative RI scheme 

In the RI scheme p6|p[ , renormalisation condi- 
tions arc imposed on amputated Green functions 
between quark states of momentum p in Landau 
gauge, setting them equal to their tree-level val- 
ues and thereby maintaining Ward identities. Eg 
for Or = (?r<7, the amputated Green function is 



A (pa) = S q (pa) 1 G {pa)S q {pa) 1 



(19) 



where the Green function, Go{po), and quark 
propagator, S q (pa), are computed by Monte 
Carlo. Then the renormalisation condition is 

Z*Va)^(HTrP A (p a )| p2=Ai2 = 1, (20) 

where Po projects onto the tree-level ampu- 
tated Green function and Z q is the wavefunction 
renormalisation constant. In this way, gauge- 
dependent non-perturbative renormalisation is 
achieved at the lattice scale /i ~ a -1 ~ 2 — 4 GeV. 

In order to obtain results in the MS scheme, 
a continuum perturbative matching calculation 
is required. This has now been done at next-to 
next-to-leading order (NNLO) 0, ie to 0(a 2 s ), 
and, when combined with the non-perturbative 
RI result, cancels the gauge dependence. 

As mentioned above, this method has been ap- 
plied to tree-level 0(a)-improved Wilson quarks 
at P = 6.0, 6.2 and 6.4, using both the vector and 
axial Ward identity definitions of quark mass j|, 
[34j. These definitions give consistent results in 
the non-perturbative RI scheme, and show no 
lattice-spacing dependence over the range studied 
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(although they are not consistent with the results 
of perturbative matching), so the best estimates, 
obtained by averaging the data, are Q 



MS 
L ud 



mf s (2 GeV) 



(2GeV) = 5.7±0.5±0.8±0.8MeV(21) 



130 ± 8 ± 15 ± 15 MeV (22) 



where the errors are due to statistics, non- 
perturbative renormalisation and systematics, re- 
spectively. Clearly, the errors in these estimates 
are too large at present to resolve the discrep- 
ancy between Eq @ (QCDSF) and Eq (|) (CP- 
PACS). 

Kilcup and Pekurovsky |3^] have reported pre- 
liminary results for the strange quark mass ob- 
tained using non-perturbative RI renormalisation 
for staggered quarks, extrapolated to the contin- 
uum from data at j3 = 6.0, 6.2 and 6.4. They 
obtain mp(2 GeV) = 129 ± 23 MeV, in good 
agreement with other estimates. 

What is needed now is a higher statistics de- 
termination, using an improved action (0(a)- 
improved Wilson, or staggered), which allows 
for a reasonably confident continuum extrapola- 
tion to compare with the CP-PACS and QCDSF 
perturbatively-matched results. 

5.3. Non-perturbative running of the 
quark mass 

As reported last year |59), the Alpha Collabo- 
ration has been implementing a non-perturbative 
method for running the quark mass in the 
Schrodinger Functional (SF) scheme from low to 
high scales. The Q(q)-improved PCAC mass in 
the SF scheme is [M 



m SF (2L) 



Z A (l + b A am q ) 



PCAC 
'lat 



(23) 



Z P (2L){1 + bp am q ) 

where w^ AC is defined by the ratio of matrix 
elements in Eq (^) , in which the axial vector cur- 
rent and pseudoscalar density are improved to or- 
der a, as in Eq (pT[). Here the scale dependence 
comes through Zp(2L), which is defined in the 
SF scheme in a box of linear size 2L. 

Non-perturbative running is implemented via 
step scaling functions, 



g 2 SF (2L) = a(2,f SF (L)) 
Z P (2L) = ap(2,g 2 SF (L))Zp(L), 



(24) 
(25) 



" P (2.r) 



0.95 




0.85 



Figure 10. Step scaling function for the quark 
mass in quenched QCD [ffOfl. 



which can be iterated to generate a sequence of 
couplings and quark masses, 



u k =g 2 SF (2 k L), 
u k+ i = a(2,u k ), 



m k = m S F(2 k L) 
m k+ i = m k /ap(2,u k ). 



(26) 



The step scaling functions are computed for dif- 
ferent lattice spacings and extrapolated to the 
continuum limit. The step scaling function for 
the quark mass is shown in Fig [w| 

Successive application of ap yields ratios 



"^SF(imax) m SF (L max /2) 



mSF(2imax) ' TOSF(2L max ) 

m SF (£max/2 8 ) 

m SF (2L max ) 



(27) 



At the smallest coupling, the SF quark mass can 
be related to the RG invariant mass, Eq (0), 
through 



— = (2b g ) exp 



dg 



dp 

bog 



(28) 



and, hence, the evolution of the running quark 
mass between this high-energy scale and the orig- 
inal low-energy scale, 2L max , is determined. This 
universal result is shown in Fig [ll[ 
At the lowest energy, Alpha obtains 



m SF {2L max ) 



0.72 fm (29) 
= 1.166 ±0.015 ± 0.008. (30) 



2L max = 1.436 r 
M 
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100 ^[GeV] 



Figure 11. Universal running of the SF quark 
mass in quenched QCD EoT] . 
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Figure 12. Renormalisation factor for the pseu- 
doscalar density in the SF scheme at the low- 
energy scale 2L max , as a function of the gauge 
coupling in quenched QCD EJ. 



So the total renormalisation factor is 

— = 1-166(17) Z A (9o)(l + b A arn q ) 
miat Z P (g , 2L max )(l + b P am q ) 

where b A — bp is small enough to be safely ne- 
glected ||(J. The matching to the lattice scheme 
is completed by computing Zp(go,2L max = 
1.436r ) for a range of (3 values, by finding 
(f3,L/a) so that LI a — 1.436 Tq /a, using known 
values of r Q /a jy]]. The results are shown in 
Fig[lJ, where the two sets of results correspond to 
using the 1- and 2-loop estimates for the bound- 
ary counterterm, ct, which is only known pertur- 
batively. All that is required now, for example, 
is the PCAC mass, (m u + m s )i a t, from a stan- 
dard simulation, to obtain the RG invariant mass 
M u + M s . Currently, this last step is missing, 
but suitable data have been generated by both 
QCDSF and UKQCD. 

6. GLUEBALLS & HYBRIDS 

Finally, I will turn to phenomenology. 
Quenched QCD predictions are playing an im- 
portant part in glueball and hybrid meson 
searches [[l2j . The challenge is to quantify mixing 
and sea-quark effects. 

6.1. The lightest glueball 

Lee and Weingarten [Q have estimated the 
mixing of the lightest discrete isosinglet states. 



Their results for the unmixed scalar quarkonium 
mass are shown in Fig fTJL together with their best 
estimate of the ++ glueball mass, from world 
data. From the relative ordering of the contin- 
uum glueball and ss meson, it is plausible that the 
/o(1710) is predominantly glue. Lee and Wein- 
garten have computed the mixing energy, and use 
it to fit the experimental data to a crude model of 
the mixing between the glueball, and the scalar ss 
and (uu + dd)/V2 (which they denote nn) states, 
with the result 

|/o(1710)) = 0.86(5)|<7)+0.30(5)|ss) 

+0.41(9)|wT) (32) 

|/o(1500)> = -0.13(5)|<7> + 0.91(4)|ss) 

-0.40(ll)|nn) (33) 

|/o(1390)) = -0.50(12)| ff ) + 0.29(9)|ss) 

+0.82(9) \nn), (34) 

which points to the /o(1710) being 74% glueball. 
They argue that, despite being largely ss, the 
opposite sign of \ss) and \nn) in |/o(1500)) sup- 
presses its decay to KK, in line with the experi- 
mental observation. 

6.2. Hybrid mesons 

The UKQCD quenched QCD result @ that 
the lightest ss exotic meson has J PC = 1 ^ and 
a mass of 2.0(2) GeV, is supported by further 
quenched results reported by MILC this year [^5| . 
SESAM has repeated the UKQCD analysis using 
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Figure 13. Continuum extrapolation of the un- 
mixed scalar ss mass, compared with the contin- 
uum limit of the scalar glueball in quenched QCD 
and the masses of the /o(1710) and the /o(1500) 
(the error bands are due to the uncertainty in 



Wilson Nf — 2 configurations at j3 = 5.6. They 
are able to perform a linear chiral extrapolation 
in the dynamical quark mass, from data at four 
quark masses, and find that the 1 h remains the 
lightest exotic, with a mass of 1.9(2) GeV, which 
is consistent with the quenched results [Q. Of 
course, exotic mesons are expected to mix with 
four-quark states. This mixing is currently under 
investigation by SESAM. 

For the wider particle physics community, pre- 
dictions for glueballs and hybrid mesons (along 
with the light quark masses) are the most inter- 
esting output of our light hadron work, and pur- 
suing a more systematic analysis of mixing and 
sea-quark effects will be an important part of our 
future programme. 

7. CONCLUSIONS & ACKNOWLEDGE- 
MENTS 

In summary, we now have precise results from 
CP-PACS which indicate that quenched QCD 
fails to reproduce the light hadron spectrum at 
the few percent level. Even so, there may be lin- 



gering doubts about the reliability of the contin- 
uum extrapolation, as there is no confirmation 
yet from improved actions. For quark masses, 
we are on the verge of resolving discrepancies be- 
tween different simulations, which are probably 
due to the unreliability of perturbative matching. 
Techniques for non-perturbative renormalisation 
are well-developed and should be mandatory from 
now on. It is time for a definitive simulation 
of quenched QCD using non-perturbative 0{a)- 
improved Wilson, or staggered quarks and non- 
perturbative matching, to establish the quenched 
light hadron spectrum and quark masses once and 
for all. 

Simulations with two flavours of dynamical 
quarks are showing hints that the strange hadron 
spectrum is in better agreement with experiment 
than in quenched QCD, although the quenched 
approximation must still be used for the strange 
quark. The biggest effect observed so far is that 
the light quark masses are around 40% smaller 
than in quenched QCD. String breaking has not 
been seen yet, but it should be within reach of 
present day simulations, and there is much to be 
done to improve our understanding of mixing and 
sea-quark effects on glueball and hybrid states. 

Somewhat surprisingly, quenched QCD has 
turned out to be a rather accurate effective the- 
ory for many quantities. Knowing its limitations 
should give us greater confidence in using it as a 
phcnomcnological model for quantities where 10- 
20% accuracy is useful. To improve on it, will 
require simulations with several different combi- 
nations of sea-quark flavours, in order to be able 
reliably to interpolate/extrapolate to the flavour 
content of the real world. We have reached the 
point where quenched QCD can be simulated 
with sufficient accuracy and this gives us a solid 
foundation from which to tackle the flavour de- 
pendence of QCD. 
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